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We investigate the tunnel rates and energies of excited states of small numbers of electrons in 
a quantum dot fabricated in a Si/SiGe heterostructure. Tunnel rates for loading and unloading 
electrons are found to be strongly energy dependent, and they vary significantly between different 
excited states. We show that this phenomenon enables charge sensing measurements of the average 
electron occupation that are analogous to Coulomb diamonds. Excited-state energies can be read 
directly from the plot, and we develop a rate model that enables a quantitative understanding of 
the relative sizes of different electron tunnel rates. 



I. INTRODUCTION 

Many readout schemes for solid-state quantum com- 
puting architectures rely on accurate control and real- 
time measurement of electron tunneling. 1 5 In GaAs 
quantum dots, tunneling rates have been tuned by con- 
trolling gate voltages as well as by exploiting energy- 
dependent tunnelingP enabling the identification of both 
orbital and Zeeman excited states. 7 9 Single-shot read- 
out using spin-to-charge conversion has been performed 
by using a capacitively cou pled qu antum point contact 
(QPC) as the charge sensorPMffi] 

In silicon, the effects of energy dependent tunneling 
are expected to be enhanced compared to GaAs, because 
of the carriers' larger effective mass. Time-averaged 
measurements of spin-dependent electron tunneling into 
quantu m dot s have been made using quantum point 
contacts! 13 * 14 * Single electron transistors have been used 
to measure tunnel rates^and to perform single-shot spin 
readout of electrons on individual dopants. 16 

Here we investigate the tunnel rates and energies of 
excited states of small numbers of electrons in a quan- 
tum dot fabricated in a Si/SiGe heterostructure. We 
find tunnel rates for loading and unloading electrons are 
strongly energy dependent: states below the Fermi level 
load more slowly as their energy decreases, and states 
above the Fermi level unload more rapidly as their en- 
ergy increases. Single-shot measurement of these tunnel 
events is achieved with good charge sensitivity. We fur- 
ther show that loading and unloading tunnel rates vary 
significantly between different excited states. This phe- 
nomenon enables charge sensing measurements to pro- 
duce plots analogous to Coulomb diamonds. The use of 
charge sensing to create such plots enables calibration 
of the gate voltage-to-energy ratio a without the need 
for measurable transport through the quantum dot it- 
self. Excited-state energies can be read directly from the 
plot, and we develop a rate model to extract quantitative 
relations between the tunnel rates from the experimental 
measurements. A simulated map of charge sensing mea- 
surements as a function of source-drain and gate voltages 
agrees well with the experimental data. We find signifi- 



cant variations in tunnel rates that may be useful in the 
loading and measurement of spins in quantum dots. 



II. METHODS 

The sample used in this work was grown by chemical 
vapor deposition on a Si(OOl) substrate with phosphorus 
doping of 1-10 Q • cm, which was polished 2° towards 
[010]. 17 The structure was step-graded to the compo- 
sition Sio.64Geo.36- The Ge concentration was then re- 
duced to Sio.68Geo.32 and a 1 fim buffer layer was grown, 
resulting in a final strain relaxation of 95%. On top of 
the relaxed SiGe, we grew 18 nm of Si (the quantum 
well), 22 nm of undoped Sio.6sGeo.32, 1 nm of doped 
Sio.68Geo.32, 45 nm of undoped Sio.6sGeo.32, and a 9 nm 
Si cap layer. The top gates were formed by electron 
beam evaporation of Pd onto the HF-etched surface of 
the heterostructure.^^ The gates sit on a square mesa 
of width 35 jam that was defined by reactive ion etch- 
ing. For the experiment reported here, the sample was 
illuminated with red light for 10 s at a temperature of 
4.2 K, before cooling to a refrigerator base temperature 
of 20 mK. Magnetoresistance measurements obtained 
from the same heterostructure give a carrier density of 
5.15 x 10 11 cm -2 and a mobility of 120, 000 cm 2 /Vs (after 
illumination). Data are reported in t his p aper with the 
quantum dot in both the one-electro n 17 l 2Q l 21 l and many 
electron regimes. 



III. RESULTS 

Fig. [TJa) shows a scanning electron micrograph of the 
gate structure of a device identical to the one reported 
here. The quantum dot in this experiment was formed 
by applying negative voltages to gates L, M, R, and T. A 
charge sensing QPC was formed by applying a negative 
voltage to gate Q L ^^The QPC was biased at 100 /iV, 
and the drain current was monitored with a low-noise 
current preamplifier. For the data shown in Figs. [T] and 
[2j the quantum dot barriers and gate voltages were tuned 
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FIG. 1. (Color online) (a) Scanning electron micrograph 
of the quantum dot top-gates, with gate labels shown. For 
pulsed-gate measurements, a pulse sequence was applied to 
gate R. (b) Chemical potential diagrams for tunneling be- 
tween the left reservoir and the quantum dot. Loading (top 
diagram) and unloading (bottom diagram) are achieved by 
varying the gate voltages Vm or Vr. (c) Typical oscillo- 
scope traces of /qpc versus time, showing real-time tunneling 
events. Downward edges indicate loading, upward edges in- 
dicate unloading. The three traces correspond to the three 
energy configurations shown in (b). (d) Fractional dot occu- 
pation as a function of Vm , obtained from traces like those in 
(c). 



such that tunneling occurred predominately through the 
left barrier, and the dot occupation could be tuned be- 
tween zero and one. 2qpc was high when the dot was 
empty and low when the dot was occupied. 

As shown in Fig.[TJc), monitoring the QPC current as 
a function of time enables measurement of single-electron 
tunneling on and off the dot. In the absence of a source- 
drain voltage Vsd across the dot, the charge occupation is 
determined solely by the alignment of the electron chem- 
ical potential in the quantum dot, which is controlled by 
the voltage Vm on gate M, with the Fermi level of the 
leads, as depicted in Fig.PRb). Three typical time traces 
of Iqpc are shown in Fig.mc), corresponding to the three 
energy configurations in panel (b). In the top trace of 
panel (c), the chemical potential of the dot is slightly 
above the Fermi level, and the dot is almost always un- 
loaded, except for occasional thermal fluctuations. The 
bottom trace corresponds to the opposite situation, and 
the middle trace corresponds to the case where the chem- 
ical potentials of dot and lead are nearly aligned, leading 
to an average dot occupation of 50%. Based on a signal- 
to- noise ratio of 32 and a preamplifier rise time of 3 ms, 
we estimate a charge sensitivity of 2.9 x 10 -3 e/y/Hz for 
our QPC measurements. 

As shown in Fig. [ljd) , data of this type enables a 
measurement of the Fermi-Dirac distribution for the dot 
occupation (ti(Vm)) = /d(m^), where ji g is the chem- 
ical potential of the dot ground state. By repeating 
measurements like that shown in panel (d) for a se- 
ries of increasing temperatures, we acquire a set of data 
that can be used to determine the proportionality con- 
stant a between the chemical potential and Vm for the 
one-electron dot!^ From a global fit we determine a = 
0.1 29 ±0.004 meV/mV, where the error bar is determined 
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FIG. 2. (Color) Pulsed gate voltage measurements, (a) A 
typical measurement of /qpc (red curve), together with its 
contribution from cross talk (dashed black curve). The volt- 
age pulse sequence applied to gate R is shown as the solid 
black curve (arb. units). Insets: chemical potential diagrams 
for loading and unloading, (b) Cumulative counting statistics 
for loading events vs. time with V pp = 8 mV. (c) Average 
of 600 current traces (blue curve) and an exponential fit (red 
curve), (d) Average of 512 current traces (red & blue curves), 
for two different values of V pp (the corresponding pulse se- 
quence is shown in black with arb. units), (e) Loading time 
vs. pulse amplitude, showing a strong increase in tunneling 
time as the quantum dot chemical potential is decreased. 



by the quality of the fit. The electron temperature is es- 
timated as 298 mK. 

By applying pulsed gate voltages to gate R, we can 
repeatably load and unload the quantum dot, enabling 
measurement of the loading and unloading tunnel rates 
as a function of the quantum dot chemical potential. A 
typical pulse sequence is shown in Fig. [2|a) (solid black 
trace). Loading occurs during the high portion of the 
gate voltage cycle, while unloading occurs during the low 
portion, as indicated in the insets. The tunnel barriers 
can be controlled independently of the dot chemical po- 
tential, although cross talk is always present. Here, we 
tune the tunnel barriers by adjusting the voltage, pri- 
marily on gate L, so that the tunneling times are much 
longer than the rise time of the current preamplifier in 
the QPC circuit. After this initial adjustment, gate L is 
held fixed, while the voltage is varied on gate R. 

A typical measurement trace corresponding to loading 
and unloading the dot is shown in red in Fig.[2|a). Steps 
in the traces occur in pairs; the first step edge is precisely 
correlated with the voltage pulse and is the cross talk 
response of the QPC to the voltage pulse on gate R. The 
cross talk component of /qpc is sketched as a dashed 



black line in Fig. [2^a) . The second step edge occurs in 
between the voltage pulses. This signal corresponds to a 
tunneling event — either the loading or unloading of the 
dot. By recording the time interval between the voltage 
pulse and the charging event, we obtain a direct, single- 
shot measurement of the loading time. 

Figure [2^b) shows a plot of the cumulative number of 
loading events as a function of the loading time, obtained 
from a pulse sequence of identical load-unload cycles with 
Vp P = 8 mV. Loading times are extracted from /qpc(£) 
in a two step process: first, we check that the dot is 
unloaded immediately prior to the loading pulse; second, 
we extract the first time that /qpc returns to the level 
corresponding to a loaded electron. The results shown in 
Fig. |5|b) take into account 584 out of 600 such current 
traces; 16 traces did not have current levels that could 
be mapped to the analysis pattern and were ignored. As 
expected, the counts saturate exponentially, with a rise 
time of tl = 82 ± 9 ms. Figure [2jc) shows a direct 
average of all 600 current traces. The decay is well fit 
by an exponential with time constant tl = 94 ± 7 ms, 
comparable to the value obtained by event counting. 

We measure the loading and unloading rates over a 
wide range of V pp . We use a pulse sequence with an 
asymmetric duty cycle, shown in black in Fig. |2|d), be- 
cause the unloading times are significantly faster than 
the loading times. The red and blue curves in Fig. [2^d) 
correspond to two different values of the pulse height Vp P 
and are averages of 512 individual /qpc traces. Loading 
and unloading times are obtained by performing expo- 
nential fits to the averaged data. Fig. [2|e) shows the 
resulting loading times as a function of Vp P . We observe 
an approximate five-fold increase in the loading time as 
Vp P varies from 1 to 6.5 mV, corresponding to a shift of 
about 343 fieV in the dot chemical potential during the 
loading phase of the cycle. 

The dependence of the tunnel rate on gate voltage is 
consistent with a simple model of transpor t in a qu antum 
dot known as energy-dependent tunneling! 6 * 9 * 11 * 25 ^ In this 
model, the tunnel rate depends exponentially on the rel- 
ative height of the tunnel barrier compared to the chem- 
ical potential of the dot, with a larger chemical potential 
corresponding to a higher tunnel rate. This explains the 
trends of the tunnel rate with respect to variations in the 
gate voltage. It also explains why the unloading rates are 
much faster than the loading rates in our dot. It is im- 
portant to note that cross talk between the plunger gate 
and the tunnel barrier would tend to have the opposite 
effect as energy-dependent tunneling, by raising the tun- 
nel barrier and lowering the tunnel rate. 

Up to this point, we have studied the tunnel rate be- 
tween the lead and the ground state of the dot, which 
varies as a smooth function of the plunger gate volt- 
age, as indicated in Fig. [5|e). We now describe how 
the tunnel rates may change more abruptly, when ex- 
cited states enter the bias window. These excited states 
may have tunneling matrix elements much larger than the 
ground state, which depend on the shape and the sym- 
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FIG. 3. (Color online) (a) Coulomb diamond plot of current 
/dot through the quantum dot in the many electron regime, 
(b) Schematic representation of two Coulomb diamonds, (c) 
Energy level diagrams showing the dot chemical potential and 
the Fermi level of the leads for two of the edges in (b). (d, e) 
Expected transconductance (solid lines) of the charge sensing 
QPC when the D and S tunnel barriers, respectively, are the 
transport bottleneck, (f, g) Color scale plot of the transcon- 
ductance g of the QPC obtained as a function of Vm and 
Vsd- Peaks in g represent transitions where dot occupation 
changes. Vr = —0.42 and —0.44 V in (f ) and (g) respectively. 
Vm differs in the two panels, to compensate for the difference 
in Vr. 

metry of the wavefunctionsJ^^^In turn, this opens up 
new doors for measuring the tunnel rates. Specifically, 
it allows us to replace gate pulsing with time-averaged 
current measurements. Time-averaging is normally as- 
sociated with methods where the current flows directly 
through the quantum dot. One example is the so-called 
Coulomb diamond plot, where the transport current is 
measured as a function of source-drain bias Vsd and the 
plunger gate voltage.^ Such transport measurements are 
challenging in the few-electron regime, where charge sens- 
ing is the preferred experimental technique!^ We show 
here that charge sensing techniques enable spectroscopy 
of silicon quantum dots, with direct correspondence to 
Coulomb diamond measurements. Similar measurements 
were previously obtained in GaAs by Schleser et alW 

Figure [3^ a) shows a Coulomb diamond plot of the cur- 
rent through the dot as a function of Vsd and Vl, for the 
case of many-electron occupation and large current flow 
through the quantum dot. Current flows only when the 
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chemical potential of the dot /i g , which depends linearly 
on Vl, lies between the Fermi energies of the source and 
drain (/is and /ip, respectively). The black regions in 
Fig. [3^ a) are blockaded. As /i g is lowered, it becomes 
level with /is; this condition corresponds to the edge of 
the diamond with positive slope (blue line in Fig. [3jb)), 
and current begins to flow. When /ijj < /i g < /i^, current 
flows, and the dot occupation alternates between TV and 
TV — 1. The condition /i g = /id corresponds to the edge of 
the diamond with negative slope (red line in Fig. |3jb)). 
When /i g < the electron does not have enough en- 
ergy to exit the dot. In this case, the electron occupation 
becomes fixed at value of TV, and the current is blocked. 

Because of the variation in the tunnel rate with gate 
voltage and the available energy levels, the same informa- 
tion available in a Coulomb diamond transport measure- 
ment can be derived from time-averaged charge sensing 
measurements. Here, we study this effect in the several- 
electron regime. Figures J3^f) and (g) show plots of the 
transconductance g = <9/qpc/<9Vm- When monitoring 
the charge sensing QPC current in this way, the physi- 
cal picture of Coulomb blockade is unaffected; however, 
instead of measuring the current through the quantum 
dot, which is immeasurably small in this regime, the 
QPC measures the average dot occupation. In the regime 
where dot occupation alternates between TV — 1 and TV, 
the average occupation depends on the electron dwell 
time, which depends on the loading and unloading tunnel 
rates. 

The interesting case for charge sensing is when Vsd ^ 
and the dot chemical potential is in the bias window. 
In this regime, the average fractional occupation /, the 
fractional part of the dot occupation (n), is given by 
/ = Tgs/(T g s + Y g c>), where Y g s is the tunnel rate from 
the source to the TV-electron ground state, and Y g r> is the 
corresponding rate to the drain. At the location of the 
purple stars in Figs. |3jd)-(g), if the tunnel barriers are 
tuned such that Y g s — Y g £>, the average fractional occu- 
pation will be / ~ 0.5. However, it is more common that 
the tunnel rates are dissimilar, and when one barrier is 
much larger than the other, it becomes a bottleneck for 
transport. 

For the limiting case T g s ^> Y 9 d, shown in Figs. [3ld) 
and (f), the electron tunnels into the dot very quickly 
when /i g < /is, but it tunnels out very slowly, so the 
average occupation will be ~ TV. A filling transition 
is therefore observed when /i g = /is, but not when 
= Vd- Thus, the transition occurs along the edge 
of the Coulomb diamond with positive slope. In the op- 
posite limit, T g D ^> Y g s, shown in Figs. 3^e) and (g), and 
achieved by changing gate voltage Vr, the dot empties so 
quickly that the average occupation is ~ (TV — 1). In this 
case, the filling transition occurs when the dot chemical 
potential is aligned with the drain, corresponding to the 
edge of the Coulomb diamond with negative slope. As 
is clear from Figs. J3^f) and (g), the tunnel barriers are 
easily tuned into either regime. 

Positive and negative slopes can also be observed with- 
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FIG. 4. (Color online) Excited-state spectroscopy, (a) Color- 
scale plot of g as a function of Vsd and Vm- The pattern in 
the transconductance is analogous to a traditional Coulomb 
diamond plot of the differential conductance through a quan- 
tum dot. (b) Schematic representation of the "transition 
map" observed in (a), showing the role of excited states and 
how the energies of these states can be extracted in analogy 
with a Coulomb diamond plot, (c) Theoretical reconstruc- 
tion df /8Vm of the transition map, based on fitting to a rate 
equation model described in the Appendix, demonstrating the 
robustness of the interpretation of the data, (d) Schematic en- 
ergy level diagrams of the dot at points A, B, and C in panel 

(b). 



out retuning the tunnel barriers, simply by expanding 
the bias window, as shown in Fig. 4^a). Here, we have 
chosen the tuning Y g s > Y 9 d. The mapping between 
the filling transitions and the conventional Coulomb dia- 
monds is sketched in Fig.Qb), where one of the diamonds 
has been shaded gray. We observe sharp corners in the 
transition map, which we attribute to excited states en- 
tering the bias window, similar to lines corresponding to 
excited states in a Coulomb diamond plot. Three such 
transitions are labelled A, B, and C in panel (b). Near 
^sd = 0, the charge transition line has a positive slope, 
since the D barrier forms the bottleneck, and the fill- 
ing transition occurs when the dot level is aligned with 
the source. In this range of Vsd, all excited states lie 
outside the bias window, and the average fractional oc- 
cupation / = T g ]j/(Tg]j + T g s) — 0. (See Appendix 
for details.) When an excited state x enters the win- 
dow, such as at point A in Fig. |4jb), the occupation is 
given by / = (T gD + Y xD )/(Y gD + Y gS + Y xD ) ~ 1. In 
this case, the approximate equality holds because of the 
strong tunnel coupling between the drain and the excited 
state: Y x r> ^> T^s ^ Y 9 d- Thus, a switch in occupation 
is caused by a switch in the rate-limiting tunnel barriers. 
We can understand the changes in slope near points B 
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TABLE I. Excitation energies and tunnel rates as extracted from Fig. [4] Here, ij specifies the tunnel channel, with i 
corresponding to the quantum dot orbital, and j corresponding to the source S or drain D lead. The tunnel rates are obtained 
by fitting the transition map in Fig.[4ja) to a finite temperature rate model described in the Appendix. The excitation energies 
(1-5) are measured relative to the ground state (g). This procedure provides tunnel rates that are normalized to the gS tunnel 
channel. The tunnel rates are then calibrated with respect to the gD tunnel channel, as described in the text. 
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and C by similar arguments. The transition map expands 
in piece-wise fashion, with sharp changes in slope indi- 
cating the presence of excited states that change which 
barrier corresponds to the bottleneck tunnel rate. Such 
changes in slope are visible over a wide range in Vsd and 
Vm- Excited states with tunnel couplings that do not 
alter the time averaged dot occupation will be invisible. 

We can analyze the experimental data of Fig. |4^a) 
using a finite temperature rate model, as described in 
detail in Appendix A. In brief, the rate equations are 
obtained by first expressing the currents for individual 
transport channels. For example, the tunnel current from 
the source to the ground state of the dot is given by 
I gS = -eT g sfs(Vg)(l - /), where /s(/x g ) is the Fermi 
function for the source reservoir, evaluated at the chem- 
ical potential of the dot ground state. By enforcing cur- 
rent conservation and assuming an infinite decay rate for 
excited states, we obtain an equation for / as a function 
of tunnel rates, excited-state energies, and temperature. 
It is important to note that, although an electron from 
the drain can tunnel into an Af-electron excited state, it 
cannot tunnel out through the same channel. Instead, 
the electron decays very quickly and tunnels out through 
the ground stated The decay process causes an intrinsic 
asymmetry, one which favors loading rather than unload- 
ing. 

The quantity df/dVyi corresponds directly to the 
transconductance of Fig.Qa). We fit this result to the ex- 
perimental data, obtaining the parameters shown in the 
next-to-last column of Table I, as well as the theoretical 
reconstruction of the transition map shown in Fig. [4jc). 
In the fitting procedure, no correlations were assumed be- 
tween the excited states in the different transport chan- 
nels. We estimate an uncertainty of about 60 /ieV for the 
excited-state energies, and an uncertainty of about 20% 
for these relative tunnel rates. 

Since transport currents are immeasurably small in 
this experimental configuration, we are not able to deter- 
mine the tunnel rates directly. In the next-to-last column 



of Table I, the tunnel rates are normalized relative to the 
gS tunnel channel, as described in the Appendix. Here, 
the tunnel rates for the N ^ (N — 1) and N ^ (N + 1) 
processes are normalized independently, since they in- 
volve different gate voltages, and should not be exactly 
equal. 

We were able to obtain an approximate calibration for 
the tunnel rates, however, by retuning the quantum dot 
to allow pulsing experiments. Shortly after the data in 
Fig. [4] was obtained, the right tunnel barrier (R) was 
pinched off, allowing electrons to tunnel through just the 
left (D) barrier. The tunnel barrier gate voltage Vl was 
left essentially unchanged. Tunnel rates corresponding 
to the gD process were then acquired on two successive 
days, with variations on the order of 25%. In this way we 
obtain the estimates T gD = 5.4 xlO 5 s _1 and8.5xl0 5 s _1 
for the JV f> (N — 1) loading and unloading processes, 
respectively. The difference between loading and unload- 
ing rates is comparable to those described above, and it 
is also consistent with predictions of energy dependent 
tunneling. In the final column of Table [TJ we use the 
mean tunnel rate estimate 7.0 x 10 5 s _1 to calibrate the 
relative rates listed in the next-to-last column. The same 
calibration provides a rough estimate for the N (iV+1) 
tunneling process, and we provide the corresponding re- 
sults in Table U 

The lever arm a, which converts gate voltage to dot en- 
ergy, is usually extracted from a Coulomb diamond plot. 
Here, we can extract a = 0.095 ± 0.004 meV/mV from 
Fig. [4^a), even though transport through the quantum 
dot itself is immeasurably slow. The uncertainty in a 
is determined by the resolution of the data plot in both 
Vsd and Vm • This value of a is different than that quoted 
above when the quantum dot was occupied by a single 
electron, because the dot has been retuned and now is 
occupied by several electrons. 
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IV. SUMMARY 



(a) 



(b) 



In summary, tunnel rates in and out of a few-electron Si 
quantum dot were measured by single-shot charge sens- 
ing. The rates were shown to depend strongly on gate 
voltage, in a manner opposite to that expected from cross 
talk, and consistent with energy-dependent tunneling. 
We also have shown that a map of the time averaged 
dot occupation, obtained by charge sensing, provides di- 
rect spectroscopic information about the quantum dot. 
Further, energy calibration — the determination of a 
— does not require retuning of the tunnel barriers, be- 
cause both the positive and negative slopes, which are 
usually extracted from Coulomb diamonds, are visible in 
a plot like that shown in Fig. [4^a). The sharp corners 
observed in the transition map are attributed to excited 
states entering the bias window and their effect on the 
tunneling bottleneck. Sudden switching of the bottle- 
neck is expected, because the tunneling matrix elements 
for excited states can be large, and because there is a 
built-in asymmetry in the tunneling process, due to the 
fast decay of excited states. 
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Appendix A: Rate equations for the transition map 

In this Appendix we present a more detailed explana- 
tion of the Coulomb diamond-like behavior observed in 
the transconductance data of Fig. (4^a). As explained in 
the main text, the presence of switches in the transition 
map indicates that the source (S) and drain (D) tunnel 
barriers are asymmetric. A change of slope occurs when- 
ever the slow barrier (the bottleneck) switches from S 
to D, or visa versa, due to excited states entering the 
bias window P Here, we will derive a simple theoretical 
description of the switching observed in this transition 
map. 
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FIG. 5. (Color online) (a) Ground (g) and excited (x,y) state 
transitions in a quantum dot, between the electron occupa- 
tions (TV — 1) and TV. The chemical potential describing tran- 
sitions between the ground states is \i g . Transitions to excited 
states may occur in either direction: when the dot loads {ji x ) 
or unloads an electron (fly), (b) The types of tunneling pro- 
cesses considered in this work. Fast decay of excited states 
essentially prohibits x-type unloading or y-type loading pro- 
cesses. 



The transition between (TV — 1) and TV-electron states 
in a quantum dot may involve excited orbitals in either 
the (TV— 1) or TV-electron manifolds. The chemical poten- 
tials for several such processes are sketched in Fig. [5ja). 
There are two necessary ingredients for observing a single 
change of slope. (1) One of the barriers (say, D) must 
form a bottleneck for transitions between the (TV — 1) and 
TV-electron ground states; in other words, T 9 d <C T^s, 
where V signifies a tunnel rate. (2) There must be an 
excited state (not necessarily the lowest excited state) 
for which the tunnel rate to D is faster than the ground 
state tunnel rate to S\ in other words, ^> Y g s. The 
fact that the decay rate T xg between the excited and 
ground states of the quantum dot is much faster than any 
of the tunnel rates also facilitates switching. Additional 
changes of slope may occur when new excited states enter 
the bias window. Note that energy-dependent tunneling 
is not a leading order effect in the switching behavior ob- 
served in Fig. [4ja) and it will not be explicitly considered 
here. 

We will consider each of the different, sequential tran- 
sition processes shown in Fig. |5|b). Although the case 
Us > I^d is presented here, analogous processes are also 
present when jis < I^d- We will adopt the notation x to 
refer processes between the (TV — 1) -electron ground state 
and an TV-electron excited state, as sketched in Fig.[5|a). 
Only one such transition is shown in the figure, although 
there may be many, in practice. The notation y refers 
to processes between the TV-electron ground state and an 
(TV — l)-electron excited state. The transition between 
ground states, g, can play a role in both loading and 
unloading of the quantum dot. However, x or y-type 
processes are essentially uni-directional, as indicated in 
panel (b), due to the fact that T xg is very large. Thus, 
when an x-type process occurs, the loaded state immedi- 
ately decays to the TV-electron ground state, before the 
dot has a chance to unload. Similarly, a y-type process 
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may unload the dot into an (TV — l)-electron excited state, 
which then decays to the ground state before loading can 
occur. Transitions between N and (N — l)-electron ex- 
cited states are not forbidden. However, they are strongly 
suppressed by the same fast decay process. 

We now derive an equation for the steady-state frac- 
tional filling of the quantum dot, /. For simplicity, 
we will assume that the filling transitions observed in 
Fig. [4ja) are thermally broadened, although lifetime 
broadening may also play a roleP^ The elastic loading 
of the quantum dot ground state from the source is then 
described by the tunneling current Is g , given by 



Here, Ts g is the tunnel rate from the source to the TV- 
electron ground state, and fs(l^g) is the Fermi function 
for the source, evaluated at the chemical potential for 
process g. The reverse process is written as 



I gS = -eT gS [l-f s (» g )}f. 



(A2) 



ISg = eT S gfs(Vg)(l ~ /)• 



(Al) 



From here on, we will assume that = T g s. 

We can write down the currents for all such processes, 
keeping in mind that a process like S —> x is allowed, 
while the reverse is not, due to the fast relaxation of the 
excited states. We can obtain a steady-state solution by 
summing up the currents through the S and D barriers 
and equating the results. Solving for /, we obtain 



TgDfD(Hg) + ^gsfs(l^g) + Ei \^ Xi D f D^Xi) + iSfs fax i )] 



FgD + T gS + Ei FxiDlDfaxi) + TxiSfsfaxi)] + E; \ T yiD (1 Id^)) + T Vi S (1 /s(/%))] ' 



(A3) 



Here, we have included all possible excited states, al- 
though some of the processes may be invisible, due to 
weak coupling to the leads. 

As an example, we consider the different cases shown 
in Fig. [6] We have focused on the bias /is > /id\ however, 
analogous results are obtained for the opposite bias. We 
will consider T g jj <C ^ g s, so that the drain barrier forms 
the bottleneck. To simplify the discussion, we will set 
T = 0, although the final fitting, shown in Fig. [4^c), 
includes temperature as a fitting parameter. 

In Figs.[6|a) and (b), we assume that /is is just slightly 
larger than /ijj. In particular, (/is — /id) < (l^g ~ Vy), 
so the g and y processes may not be in the bias window 
simultaneously. At T = 0, the arrangement shown in 
Fig. 6ja) involves no loading processes, as consistent with 
Fig. 5Tb) and the discussion in the main text. The Fermi 
functions take the values or 1, and Eq. (A3) reduces to 
/ = 0, indicating an empty dot. On the other hand, for 
the arrangement shown in Fig. |6[b), the g process can 



load and unload. In this case, Eq. (A3) reduces to 



/ = 



gs 



1. 



gS 



gD 



(A4) 



In this regime, the dot is filled because of the bottleneck 
at the drain barrier. Thus, the filling transition occurs 
when n g = /is- 

In Fig.[6jc), the bias is increased such that (/is — Hp) > 
(/ig — /i y ). As before, when /i g > /is, Eq. ( |A3[ ) reduces 
to / = 0. When /i g < /is, the quantum dot can load 
through channel g, but it can unload throug h ch annels g 
or 2/, as consistent with Fig.[5jb). Now Eq. (A3) reduces 
to 



gs 



gs 



gD 



yD 



(A5) 



where the final equality holds when the tunnel rate for 
the process y — » D is large; that is, when T y jj ^> T^. 



When /iy < /id, we recover Eq. (A4). Thus, for the ar- 
rangement shown in Fig.|6|c), the filling transition occurs 
when /i y = /ijj. The alignment of the filling transition 
with the drain rather than the source causes a change 
of slope in the transition map, as discussed in the main 
text. The switching of the slope occurs precisely when 
(/is — /id) = (l^g — l^y), corresponding to the excitation 
energy of the (N — l)-electron dot. 

In Fig.|6jd), the bias is increased such that (/is—^d) > 
(/i x — /iy). When all three chemical potentials lie inside 
the bias window, Eq. (A3) reduces to 



gS 



1. 



gs 



gD 



yD 



■ xS 



(A6) 



The final equality holds if we assume that T x s ^> T^. 
Thus, the filling transition occurs when the x process 
aligns with the source, again causing a change of slope in 
the transition map. In this case, the change of slope oc- 
curs precisely when (/is — ^d) = (Hx—^y), corresponding 
to the point where the bias is equal to the sum of exci- 
tation energies for the (TV — 1) and A/"-electron dots. 

The transition map therefore gives a direct method 
for performing spectroscopy in a quantum dot, since the 
changes of slope correspond to excited state transitions 
entering the bias window. The corresponding energies 
can be read off directly from the Vsd axis in Fig. [4] 
However, the excited states can only be observed if they 
couple strongly to the lead. 

It is not surprising that consecutive transitions should 
have larger tunnel rates, since excited states tend to be 
more spatially extended. However, it is also possible for 
excited states to be poorly coupled to the leads. Such 
states will not cause switching, and will be rendered in- 
visible in the transition map. Similarly, a given level may 
couple differently to the source and drain reservoirs. For 
one reservoir, it may form an important tunnel channel, 
while for the other, it may form a bottleneck. 
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(a) 



(b) 



(c) 



(d) 





_y_ 






















_y_ 








D 




s 




D 




s 




D 



FIG. 6. (Color online) Cases considered for discussing the quantum dot filling, Eq. (A3) 



The fitting result shown in Fig. |4^c) was obtained by 
simultaneously fitting the tunneling and energy param- 
eters, appearing in the derivative df /8Vm of Eq. (A3), 
to the transconductance data of panel (a). We have in- 
cluded a "noise floor," below which all the data were 
assumed to be indistinguishable from zero, as was done 
in the experimental plot. We also assumed a saturation 
value for df /8Vm, corresponding to the color yellow. Fi- 
nally, we included temperature as a fitting parameter. 

Table I shows the results of our fitting analysis for one 
of the Coulomb diamonds, corresponding to the shaded 
diamond in Fig.Qb). We also obtain the temperature es- 



timate of T ~ 0.6 K. This result does not match the esti- 
mate obtained from Fig.JIJd), due to differences between 
conventional and pulsing transport techniques, and the 
errors introduced by the low and high transconductance 
cut-offs used in our fitting. 

Note that the bright spots in the transition lines with 
negative slope, at the top right of the transition map in 
Fig. [4^c), are real. They occur along extensions of lines 
associated with excited states from lower diamonds, as 
shown in panel (b). The corresponding bright spots are 
offset in panel (a), due to a charging event which shifted 
the transition line slightly to the right. 
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